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The  article  presents  torque  and  flux  ripples  minimization  based  on  adaptive 
filter.  The  adaptive  filter  coefficients  optimized  and  adapted  on  line  by  using 
particle  swarm  optimization  (PSO)  technique.  The  proposed  methodology 
applied  on  closed  loop  speed  control  based  on  direct  torque  control  (DTC) 
for  surface  mounted  permanent  magnet  synchronous  motor  (PMSM).  Fuzzy 
logic  controller  (FLC)  used  as  speed  controller  while  proportional-integral 
(PI)  controller  used  as  torque  and  flux  controllers.  Second  order  infinite 
impulse  response  (HR)  filter  is  used  for  ripple  reduction  generated  due  to 
FLC.  The  drive  system  modeled  using  Matlab/Simulink  software  in  order  to 
dynamically  evaluate  the  performance  of  the  proposed  drive  system  at 
different  operating  conditions.  The  results  provide  evidence  that  the  adaptive 
filter-FLC  integration  with  optimal  coefficients  minimizes  torque  and  flux 
ripples  with  reduction  of  total  harmonic  distortion  generated  in  the  three- 
phase  currents. 

Copyright  ©  2019  Institute  of  Advanced  Engineering  and  Science. 

All  rights  reserved. 


Corresponding  Author: 

Yasser  Ahmed, 

Taif  University,  Saudi  Arabia. 

Email:  abdelsalamyasser@yahoo.com 


1.  INTRODUCTION 

One  of  the  famous  drive  systems  is  permanent  magnet  synchronous  motor  (PMSM)  that  is  used  in 
the  applications  which  require  high-performance  such  as  robotics,  transportation  and  industry.  The  main 
advantages  of  PMSM  are  high  efficiency,  small  losses,  low  motor  size,  high  power  density,  high  acceleration 
and  deceleration.  Speed  of  PMSM  can  be  controlled  by  using  either  direct  torque  or  field  oriented  control 
(DTC  or  FOC).  The  DTC  is  preferred  because  of  its  rapid  dynamic  torque  response,  simple  structure  and 
robust  control  scheme.  However,  DTC  has  drawbacks  such  as  ripples  in  torque  and  flux,  which  cause 
considerable  amount  of  noise  and  vibrations  which  affect  the  derive  performance  [1],  [2]. 

There  are  two  operating  techniques  of  the  DTC,  space  vector  modulation  (SVM-DTC)  and 
conventional  switching  table  (ST-DTC).  In  the  SVM-DTC,  arbitrary  voltage  vectors  with  variable 
magnitudes  and  fixed  switching  frequency  are  generated  for  controlling  the  torque  and  flux  with  less  ripples 
that  the  SVM-DTC  is  more  accurate  than  the  ST-DTC.  In  order  to  control  the  torque  and  flux  based  on  the 
ST-DTC,  the  appropriate  voltage  vector  is  pinpointed  from  switching  table  of  a  certain  number  of  specified 
voltage  vectors  that  are  estimated  at  variable  switching  frequency.  Therefore,  the  SVM-DTC  is  used  in  this 
work  [3],  [4], 

The  used  methodologies  for  ripples  minimization  in  the  motor  drives  are  optimal  design  of  the 
machine  and  suitable  control  scheme.  The  design  process  should  optimize  the  stator  and  rotor  parameters  to 
reduce  the  cogging  torque  as  well  as  harmonic  contents  of  the  magnetic  flux.  Further  study  was  in  [5]  in 
order  to  reduce  the  cogging  torque  where  the  proposal  was  based  on  the  design  of  slots  widths  and  tooth 
profiles.  In  [6],  the  author  studied  the  combination  influence  of  consideration  of  the  pole  and  slot  shapes  on 
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torque  ripples  in  a  PMSM  and  the  results  clarified  that  the  ripples  maybe  increased  if  the  machine  is  not 
carefully  designed. 

Many  control  methodologies  have  been  reported  in  the  literature  for  limiting  flux  and  torque  ripples 
in  which  it  maintained  smooth  operation  of  the  drive.  The  control  methodology  adapts  the  controller 
parameters  to  optimize  the  output  voltage  vector  from  the  inverter  [7],  [8].  Toward  the  speed  control  using 
FOC  in  the  interior  PMSM  drive  system,  [9]  suggested  to  use  the  adaptive  filter  in  order  to  minimize  the 
torque  ripples  where  the  fuzzy  logic  controller  (FLC)  was  designed  and  the  filter  coefficients  selection  was 
based  on  the  trial-error  concept.  The  speed  hannonic  magnitude  is  used  as  a  feedback  control  signal  in  a 
closed-loop  current  controller  in  order  to  reduce  torque  ripples.  The  obtained  results  confirm  the  attained 
ripples  reduction  [10],  In  [11],  a  modified  switching  table  is  designed  for  the  three  voltage  vectors  in  DTC 
selection  rather  than  two  voltage  vectors.  The  corresponding  less  ripple  contents  are  ascertained  in  torque  and 
flux  with  sinusoidal  stator  current.  However,  the  durations  of  the  three  vectors  are  computed  to  reduce  the 
switching  frequency.  In  [12],  variable  amplitude  and  angle  voltage  vector  have  been  controlled  to  reduce 
ripples  based  on  DTC  for  PMSM  drive  system.  A  predictive  DTC  is  proposed  to  enhance  the  conventional 
DTC  performance  during  transient  and  steady  state  conditions.  For  fast  transient  response,  a  parameter 
regulation  of  the  voltage  vectors  is  considered.  However,  the  ripples  are  minimized  by  investigating  the 
optimal  amplitude,  angle,  and  duration  of  the  voltage  vector  [13]. 

Artificial  intelligence  (AI)  approaches  were  used  for  optimal  selection  of  control  scheme  parameters 
to  minimize  torque  and  flux  ripples  and  enhance  system  performance  as  reported  in  [14].  The  famous 
examples  of  Al-based  methods  are  genetic  algorithm  (GA)  as  well  as  particle  swarm  optimization  (PSO). 

In  this  work,  the  speed  of  surface  mounted  PMSM  controlled  using  DTC  technique.  FLC  integrated 
with  second  order  infinite  impulse  response  (HR)  filter  used  as  speed  controller  while  PI  controllers  used  in 
torque  and  flux  loops.  The  HR  filter  used  to  reduce  the  ripples  generated  by  FLC.  PSO  is  used  to  optimize  the 
system  performance  on-line  by  selecting  the  optimal  coefficient  of  adaptive-IIR  filter  and  controllers  which 
minimize  flux  and  torque  ripples.  The  adaptive  filter  coefficients  are  optimized  simultaneously  and  varied  as 
a  function  in  the  torque  ripples.  A  closed  loop  speed  control  based  on  FLC  and  SVM-DTC  in  surface 
mounted  PMSM  is  proposed.  The  proposed  system  modelled  in  Matlab  environment.  The  proposed  method 
is  justified  by  applying  the  methodology  on  different  operating  conditions  of  the  system  under  investigation. 


2.  TEST  SYSTEM  STRUCTURE 

The  presented  test  system  consists  of  PMSM  with  inverter  bridge,  speed  FLC  integrated  with  ripples 
minimization  algorithm  and  DTC  controller  as  depicted  in  Figure  1.  The  system  components  are  as  follows. 


2.1.  PMSM  mathematical  model 

The  PMSM  model  in  rotating  d-q  frame  is:  [15],  [16]: 

ud=  Rs.id+^f-wripq  (1) 

uq  =  Rs.iq+^-  +  wripd  (2) 

H^d  ~  id  4"  *Pm  (3) 

^Pq  ~  Lq-iq  (4) 

Te  =\p  ('Pd-iq-'Pq-id)  (5) 

From  3,  4  and  5,  the  electromagnetic  torque  ( Te )  can  be  derived  as: 

—  ~P  (jpm‘  iq  T  (/' f /  “  f (/ ) k/ ■  id)  (6) 


From  (6)  shows  that  Te  contains  two  components  that  are  the  torque  due  to  permanent-magnet  and 
the  reluctance  torque.  For  the  PMSM  type  considered  in  the  study,  the  surface  mounted  type  has  thatL^L^. 
This  mean,  the  reluctance  torque  is  eliminated,  Te  equation  can  be  written  as: 

Te=\p^miq  (7) 


Adaptive  filter-flc  integration  for  torque  ripples  minimization  in  pmsm  using  PSO  (Yasser  Ahmed) 


50 


ISSN:  2088-8694 


□ 


Te="P\'l>s\\'l>m\sin8  (8) 

The  mechanical  equation  of  the  PMSM  can  be  formulated  by  [17]: 

Te  =TL+Jm^+Bmwr  (9) 

where  ud,  uq,  id  and  iq  are  stator  voltages  and  currents  respectively  in  d-q  rotating  frame.  Rs  is  the  stator 
resistance.  wr  is  the  mechanical  speed,  i pd,  xpq,  Ld,  and  Lq  are  stator  flux  and  inductance  components, 
respectively.  ips  and  i/im  are  the  stator  flux  linkage  and  permanent  magnetic  flux,  respectively.  P  is  the  motor 
pole  pairs  number.  5  is  the  load  angle.  TL  is  the  load  torque.  Bm  and  ]m  are  friction  coefficient  and  composed 
motor-load  inertia,  respectively. 


2.2.  Speed  controller 

The  FLC  acts  as  a  PI  speed  controller  in  the  DTC  drive  system  as  depicted  in  Figure  1.  By 
differentiation  of  9,  the  speed  is  regulated  based  on  controlling  the  change  in  electromagnetic  torque  (AT;) 
that  is  presented  in  (10).  The  components  of  speed  controller  explained  in  the 
following  sections. 


at; 


ATL+Jm^+BmAwr 


(10) 


2.2.1.  FLC  model 

The  FLC  consists  of,  fuzzification,  fuzzy  inference,  and  defuzzification  as  depicted  in  Figure  2.  The 
speed  error  ( E )  and  its  change  (AE)  are  the  FLC  inputs  while  the  torque  change  (A  Te* )  is  the  output.  The 
inputs  and  output  relationships  of  the  FLC  can  be  written  as  follows: 


E  =  w;(n)  —  wr(ri) 

(ID 

AE  =  E(n )  -  E(n  -  1) 

(12) 

at;  =/(E,AE) 

(13) 

where  E(n  —  1)  is  a  speed  error  pre-sample,  E(ri),  w(!(n)  and  wr(n)  are  the  present  samples  of  speed  error, 
reference  speed  and  mechanical  motor  speed  respectively. 

The  fuzzification  process  maps  the  input  variables  E  and  AE  to  the  corresponding  linguistic 
variables.  The  fuzzy  inference  gets  the  control  rules  for  every  input.  The  defuzzification  converts  the  input 
rules  to  the  crisp  value  of  output  AT;. 
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Figure  2.  Fuzzy  logic  speed  controller 


The  E  and  EE  membership  functions  are  depicted  in  Figure  3(a).  The  designed  sets  of  the  input 
variables  are  as  that  PB,  PS,  ZE,  NS,  and  NB  are  positive  big,  positive  small,  zero,  negative  small,  and 
negative  big,  respectively.  The  output  variable  A  Te*  membership  functions  are  depicted  in  Figure  3(b)  where 
the  output  sets  are  that  PB,  PM,  PS,  ZE,  NS,  NM,  and  NB  are  acronym  of  positive  big,  positive  medium, 
positive  small,  zero,  negative  small,  negative  medium,  and  negative  big,  respectively.  Membership  types  and 
limits  are  obtained  by  experience  [18]. 


(a) 


(b) 


Figure  3.  Membership  of  FLC 


The  fuzzy  inference  system  uses  Mamdani-type  method  where  control  rules  of  the  fuzzy  controller 
done  by  if  and  then  logic  operators  according  to  the  normalized  values  of  Hand  AE  in  order  to  estimate  the 
set  of  ATg  as  listed  in  Table  1 .  The  centroid  method  is  used  as  the  defuzzification  method  to  compute  the  ATg 
value  to  reduce  the  value  of  E . 


Table  1 .  Fuzzy  logic  rules 


A  E 

~R 

NB 

NS 

ZE 

PS 

PB 

NB 

NB 

NB 

NB 

NM 

ZE 

NS 

NB 

NM 

NS 

ZE 

PM 

ZE 

NB 

NS 

ZE 

PS 

PB 

PS 

NM 

ZE 

PS 

PM 

PB 

PB 

ZE 

PM 

PB 

PB 

PB 

2.2.2.  Ripples  minimization  algorithm 

The  main  source  of  speed  and  torque  ripples  in  the  proposed  DTC  system  comes  from  FLC.  2nd 
order  adaptive  HR  filter  utilized  for  ripples  minimization.  This  proposed  model  is  as  illustrated  in  Figure  4 
where  it  can  be  represented  by  the  following  time  difference  equation  [19]: 

y[n ]  =  a0  x[n ]  +  %  x[n  —  1]  +  a2  x[n  —  2]  +  fqyfn  —  1]  +  b2  y[n  —  2]  (14) 

where  y  represents  the  filter  output  (filtered  AE),  x  represents  the  filter  input  (E)  and  n  represents  the 
signal  index. 
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The  filter  coefficients  are  optimized  using  PSO  algorithm  to  minimize  the  torque  and  flux  ripple 
contents.  Each  filter  coefficient  can  be  represented  by  at  =  kj  +  kf  *  AT,  where  AT  =  | AT,J  —Te  |.  The 
coefficients  kj  and  kf  are  optimized  simultaneously  by  PSO  over  all  the  operating  points. 


Figure  4.  HR  filter  model 


2.3.  Direct  torque  control 

The  voltage  vector  selection  is  the  main  key  in  the  DTC,  where  a  design  of  the  voltage  vector  is 
made  for  the  regulation  of  the  electromagnetic  torque  and  stator  flux.  From  (8),  the  torque  is  a  function  in  the 
5.  The  angle  5  between  ips  and  i/>m  remains  constant  during  steady  state  where  both  ips  and  xpm  run  at  the 
same  speed.  During  transient  state,  5  is  changed  because  of  the  rotation  of  both  ips and  ipm  at  different  speeds. 
Comparing  with  respect  to  electrical  time  constant,  the  system  mechanical  time  constant  is  slower  response, 
then  the  angle  5  can  be  adjusted  by  regulation  of  the  ips  speed  with  respect  to  the  ipm.  This  is  implemented  by 
selecting  the  suitable  voltage  vector  for  the  error  minimization  of  the  torque  and  flux  with  respect  to 
reference  values.  DTC  block  consists  of  Pi-torque  controller,  Pi-flux  controller  and  flux/ 
torque  estimator. 

2.3.1.  Torque  controllers 

The  torque  and  flux  controllers  utilize  conventional  PI  controllers.  The  input  of  the  torque  controller 
is  the  difference  of  the  actual  torque  with  respect  to  the  reference  value  in  order  to  attain  the  output  uq.  The 
torque  estimation  is  based  on  the  measured  current  and  estimated  flux: 

Te^QVaip-Vpia)  (15) 

2.3.2.  Flux  controller 

The  input  of  the  flux  controller  is  the  difference  between  actual  and  reference  fluxes  while  ud  is  the 
attained  output.  Then,  the  controllers  output  transformed  to  control  voltage  components  in  a-(S  frame  ( ua ,  Up) 
which  are  used  to  generate  the  SVPWM  pulses  that  fire  the  inverter  bridge.  The  generated  voltage  vector 
should  maintain  the  torque  and  flux  at  its  required  value  with  minimum  ripples.  The  magnitude  and  angle  of 
f's  can  be  estimated  by  the  following  equations  [20], 


%  =  J>T2a  +  V2p  (16) 

0s  =  tan_1—  (17) 

V a 

where: 

Va  =  I  ua  -  ia  dt  (18) 

Vp=fup-ip  dt  (19) 
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where  ia  and  ip  represent  the  current  components  in  a-/?  frame,  Va  and  Fp  represent  the  a-p  components  of 
the  estimated  Fs. 


3.  PSO  ALGORITHM 

PSO  as  one  of  meta-heuristic  optimization  techniques  depends  on  searching  the  optimal  solution 
within  the  search  area  based  on  the  exchange  of  experiences  among  particles  in  the  population.  The  particles 
in  the  swarm  modify  their  positions  in  the  next  iteration  based  on  individual  local  best  position  and  global 
best  position  of  swarm.  Each  particle  represents  a  solution  for  the  control  variables  of  the  entire  optimization 
problem.  The  control  variables  which  used  for  minimization  of  torque  ripples  in  PMSM  are  the  speed,  flux, 
torque  PI  gains  controllers  and  filter  coefficients.  For  the  particle  (i),  it  modifies  the  position  in  the  next 
iteration,  £+1,  (T'/c+1)  based  on  the  velocity  (Sf+1)  of  the  particles  using  the  following  equations  [21]: 

Pf+1  =pjc+  xSt+1  (20) 

Sf+1  =  zk+1S?  +  Ci(s£  -  Sf)  +  C2(s£  -  St)  (21) 

where,  x  is  adaptive  factor  which  used  to  control  the  swarm  conversion,  z  is  a  weight  factor  for 
controlling  the  velocity  of  swarm  towards  the  optimal  solution.  Cl  and  C2  are  generated  random  numbers  in 
the  range  of  0  and  2.  Sp  and  Sgi are  the  local  and  global  best  positions  at  the  previous  iterations. 

The  objective  function  is  to  minimize  electromagnetic  torque  ripples  (Jrlpple),  which  can  be 
represented  by  the  following  expression: 


Triple  =  UCnO  -  Ty  (22) 

T  =  “HiLiTCO  (23) 

where  n  is  the  number  of  samples,  f  is  the  average  value  of  electromagnetic  torque. 

The  optimal  solution  should  satisfy  the  operational  constraints.  The  considered  constraints  are  the 
stator  flux  ripples  ( Fripple )  and  mean  of  speed  error  (iv)  can  be  represented  by  the  following  expressions 
[22],  [23], 


^ ripple  =  J^I?=i(phi(i)  -  phi)2 

(24) 

Phi  =  phi(i) 

(25) 

w  =  ^ILiwCO 

(26) 

where,  phi  is  a  stator  flux  average  value. 

The  constraints  are  taking  into  account  during  optimization  by  an  adaptive  penalty  function  where  a 
modified  objective  function  fonnulated  based  on  the  normalized  original  objective  function  and  normalized 
constraints  violations  [24],  [20].  The  modified  objective  function  (Tf}pplef)  which  penalized  the  infeasible 
solutions  for  considering  the  constraints  violations  (Fiof£)  can  be  formulated  by: 


j>mod,k 
1  ripple 


i  =  J(T&)2  +  (VKD2  +  Ft  *  X  +  Ff  *  7ft 


(27) 


Where  Ff  and  Ft  are  the  number  of  feasible  and  infeasible  solutions  in  swarm  respectively.  X  is 
penalization  function  which  equal  zero  at  no  feasible  solutions  (Ff  =  0)  or  Vioft  with  no  infeasible  solutions 
(F;  =  0)  in  swarm  at  the  entire  iteration.  The  schematic  diagram  for  the  iterative  process  presented 
in  Figure  5. 


Adaptive  filter-flc  integration  for  torque  ripples  minimization  in  pmsm  using  PSO  (Yasser  Ahmed) 


54 


ISSN:  2088-8694 


□ 


Figure  5.  Schematic  diagram  for  PSO 


4.  SIMULATION  RESULTS 

Matlab/Simulink  software  is  utilized  in  order  to  evaluate  the  performance  of  the  proposed  drive 
system  presented  in  Figure  1.  The  PSO  algorithm  written  and  tested  using  Matlab  m-file.  The  PMSM  is 
driven  at  100  rad/sec  speed  while  a  disturbance  in  load  is  applied  from  50  to  100  N.m.  The  system  simulated 
with  and  without  ripples  minimization  algorithm,  where  PSO  used  to  attain  the  optimal  system  parameters  in 
each  test.  The  PSO  parameters  which  used  during  optimization  process  listed  in  Table  2.  The  PMSM 
parameters  and  the  obtained  optimal  filter  coefficients  are  as  illustrated  in  Table  3. 


Table  1.  PSO  parameters 


Variable  name 

value 

Variable  name 

value 

Population  size 

10 

Max.  iterations 

1000 

Ci 

1.5 

Zmax 

0.95 

c2 

1.5 

^min 

0.4 

X 

1 

No.  of  constrains 

2 

Table  2.  PMSM  parameters  and  optimal  filter  coefficients 


Filter  coefficients 

Motor  parameters 

aj,  al 

0.07,  0.01 

R  (ohm) 

0.0068 

a\,al 

0.66,  0.2 

Ld  (mH) 

0.482 

ah  a\ 

0.25,  0.22 

Lq  (mH) 

0.482 

blbl 

-0.51,-0.22 

P  (poles) 

4 

bh  hf 

-0.04,  -0.25 

J  (kg-m2) 
Ym(Wb) 

0.0015 

0.1413 

Figure  6  and  Figure  7  present  the  performances  of  the  drive  system  without  and  with  filter 
respectively.  Figure  6(a)  and  Figure  7(a)  show  the  speed  response  where  speed  reference  is  set  at  100  rad/sec. 
Figure  6(b)  and  Figure  7(b)  show  the  stator  flux  response  where  the  flux  reference  is  set  at  0.1413Wb. 
Figures  6(c)  and  7(c)  illustrate  torque  responses  where  the  load  torque  changed  from  50  to  100  Nm  at  t=0.5 
sec  (partially  loaded).  Figure  6(d)  and  Figure  7(d)  show  the  three-phase  current. 

Figure  6  shows  a  very  fast  response  during  starting  and  load  step  change  with  approximately  zero 
steady  state  error  where  the  flux  and  torque  performance  contain  considerable  ripples  and  total  harmonic 
distortion  (THD)  in  current. 

The  drive  system  performance  including  filter  effect  is  depicted  in  Figure  7  at  the  same  operating 
conditions  when  the  system  operates  without  filter.  The  torque  ripples  reduced  from  3.86  to  1.038  N.m  and 
the  flux  ripples  reduced  from  2.9  *10-3  Web  to  9.4*10-4  Web  when  using  the  proposed  ripple  minimization 
algorithm  which  is  within  specified  limit  of  5  *10-3  Web  during  optimization,  while  THD  in  current 
waveform  is  reduced  from  12%  to  4.8%. 
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Figure  6.  Drive  response  without  filter 


Figure  7.  Drive  response  filter 
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5.  CONCLUSION 

This  paper  presents  the  application  of  adaptive  HR  filter  integrated  with  FLC  as  a  speed  controller  to 
minimize  the  torque  and  flux  ripples  which  coming  from  FLC.  The  ripple  minimization  algorithm  used  to 
update  the  filter  coefficients  on  line  based  on  feedback  signal  from  the  torque  ripples.  The  optimal 
coefficients  of  HR  and  controller  parameters  have  been  estimated  by  using  PSO  to  minimize  the  torque  and 
flux  ripples.  The  system  behavior  with  and  without  the  proposed  ripple  reduction  method  has  been  tested  by 
simulation.  The  simulation  results  show  a  reduction  in  the  ripples  of  the  torque  and  stator  flux  has  been 
attained  by  73%  and  67%  respectively  with  respect  to  the  system  performance  without  ripple  reduction 
algorithm.  However,  the  total  harmonic  distortion  in  the  current  waveform  has  been  reduced  by  60%. 
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